Second-order polaron resonances in self assembled quantum dots 



Piotr Kaczmarkiewicz and Pawel Machnikowski 
Institute of Physics, Wroclaw University of Technology, 50-370 Wroclaw, Poland 

We theoretically study the optical properties of an InAs/GaAs quantum dot (QD) near the area 
of the second-order resonance between an electron confined in the QD and two longitudinal optical 
phonons. We present the absorption spectra of an inhomogeneously broadened QD ensemble and 
show that the minimal model needed for an accurate description of such a system needs to account 
for 3-phonon states. We study also the influence of the QD height to width ratio on the optical 
properties of the polaron system. The dependence of the width of the resonance and the position 
of the second-order resonant feature on the height to width ratio is presented. 
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I. INTRODUCTION 

Interaction between confined carrier states and quanta 
of lattice vibrations is important in the description of 
both optical and transport properties of nanostructures. 
Coupling of discrete carrier states localized in a quantum 
dot (QD) with nearly dispersionless longitudinal optical 
(LO) phonons leads to formation of correlated electron- 
phonon states, i.e. polaron states Such coupling of 

electronic and phononic states results in the appearance 
of anticrossings whenever the energy of an unperturbed 
carrier state approaches the energy of one or multiple LO 
phonons [l|, Q • Particularly interesting is the resonance 
of second-order, that is, the resonance between a purely 
electronic state and a state with two LO phonons. Its 
description is nontrivial since there are no direct matrix 
elements coupling those states. In our recent works 
we have shown that such a resonance can be accurately 
modeled using the Frohlich coupling Hamiltonian with 
standard (bulk) values of material parameters, contrary 
to earlier suggestions [1]. 

In this paper, we compare the calculated single QD 
polaron spectra, as well as intraband absorption spectra 
of inhomogeneously broadened QD ensembles obtained 
for different computational bases. We also study how 
the spectrum in the resonance area depends on the QD 
geometry (height to width ratio). 

The paper is organized as follows. First, in Sec. |lTl we 
present the model used for the description of an electron 
interacting with a phonon subsystem. Then, in Sec. IIII[ 
numerical results are presented. The effect of reducing 
the computational basis on the spectrum of the polaron 
is discussed, and the influence of the QD shape on the 
resonance parameters (the position and width of the res- 
onance) is studied. The conclusions are presented in Sec. 
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II. MODEL 

We consider a symmetric, self-assembled InAs QD oc- 
cupied by a single electron. The carrier confined in the 
parabolic potential of the QD interacts with the polariza- 



tion field associated with LO phonons in the surrounding 
material (GaAs). The parabolic potential model remains 
qualitatively valid as long as the confinement supports a 
sufficient number of electronic shells (3 in our model, see 
below). Otherwise, a more general model, including the 
continuum states, would need to be used. 

The Hamiltonian of the system can be expressed as 

H = Hq + Hint + Hph, (1) 

where Hq is the Hamiltonian of an isotropic QD, iJph is 
the free phonon term and Hint describes electron-phonon 
interaction. The dynamics of an electron in the direction 
perpendicular to the QD plane is restricted to the ground 
state of a harmonic oscillator with the out-of-plane con- 
finement length Iz = \J'hl {m*uJz), where m* = 0.067mc 
is the effecitve mass of the electron and hujz is the energy 
separation between eigenenergies of quantum harmonic 
oscillator in the z direction. This is a very good approxi- 
mation since the confinement in that direction is usually 
much stronger than the confinement in the xy plane, that 
is, ujQ <^ ujz, where Hujq is the energy separation between 
the in-plane states. 

The eigenfunctions of an electron in the two- 
dimensional parabolic potential are described by the 
Fock-Darwin states Q- The Hamiltonian describing the 
2-dimensional states in an isotropic quantum dot. 



can be expressed in the basis of Fock-Darwin states as 
Hq = y^e„m|nTO)(nm|, 

nm 

where n, m are radial and angular quantum numbers and 
the respective eigenvalues are 

e„m = hujQ{2n + \m\ + 1). 

An important fact is that only the LO phonons from 
the vicinity of F point are coupled to the electronic sub- 
system [8j . Moreover LO phonons are nearly dispersion- 
less for small wave vectors and it is possible to represent 
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the phonons as a set of discrete collective modes rather 
than a continuum of states ill. 

The coupling between the confined electrons and LO 
phonons in the basis of Fock-Darwin states written in 
terms of collective phonon modes has the form [5| . 
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nmn'm' a 
^"/(nm)(n'm')a-E^m' —m,a + H.C., 

where 17 is the fr equency of LO phonons at the F point, 
Iq = fi / (m* ujq) is the in-plane characteristic length, eo 
is the vacuum permittivity, e = (I/Eoo ^ I/es)^^ = 70.3 
is the effective dielectric constant, 7(nm)(n'm')a ^-re the 
coupling constants for the collective LO-phonon modes 
0, and B is collective phonon mode anihilation operator. 
The last part of Eq. (P), 

i 

describes free LO phonon term, where i numbers all dif- 
ferent collective phonon modes included in the model. 

After diagonalization of the Hamiltonian ([T]) one can 
easily calculate the intraband transition amplitudes be- 
tween polaron ground state (s-like) and higher excited 
polaron states, which in the presented range of energies 
consist mostly of polaron p-like states. 

As the intrinsic line width observed in the experiment 
[H is much smaller than the ensemble inhomogeneity as- 
sumed in our numerical computations we do not include 
line broadening effects in our model. Formally, our mod- 
eling corresponds to the zero temperature situation. 
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FIG. 1: Comparison between the polaron spectra for a single 
QD in the area of the second-order resonance obtained by 
diagonalization of the Hamiltonian for different truncations 
of basis states. The QD height to width ratio is A = 0.2. 
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FIG. 2: Comparison of absorption spectra for different trunca- 
tion of computational basis. Feature A - the main absorption 
peak, B - the two-phonon resonance, C - phonon replica of 
main absorption peak, D - phonon replica of the second-order 
resonance. The QD height to width ratio A = 0, mean transi- 
tion energy e = 63 meV, and standard deviation a — 5 meV. 



III. NUMERICAL RESULTS 

In our numerical modeling, we choose 3 lowest lying 
electronic shells (6 states). For such a number of elec- 
tronic states, 14 different collective phonon modes need 
to be constructed. We restrict the computational basis 
to states with up to 4-phonons and show that a minimal 
model describing the second-order resonant polarons has 
to include states with up to 3-phonons. Although the 
presence of the 4-phonon states has a minor impact on 
the two-phonon resonance it can be important when con- 
sidering higher order effects. 

One of the methods of formally bringing electronic and 
phononic states to resonance is QD size adjustment. The 
energy separation between the s and p shell of confined 
electron states, huJo, decreases with growing dot size and 
for a certain size the resonance condition between the s-p 
splitting and twice the LO phonon energy is met. Since it 
is the energy level splitting that is relevant for the spec- 
tral properties (rather than the QD size), it is convenient 
to use this quantity as an independent parameter in the 
model. 



In Fig. [U we present the energy spectrum for the case 
of a QD with a realistic shape (A — Iz/la — 0.2) for dif- 
ferent numbers of phonon modes taken into account in 
the model. As can be seen, there is a big difference be- 
tween the results obtained from models including 2- and 
3-phonon states, whereas adding 4-phonon states does 
not produce any important changes in the presented area 
of the polaron spectrum. If only states with up to two 
phonons are included an additional line appears in the 
area of the second-order resonance and the whole spec- 
trum is shifted towards higher energies. When 3-phonon 
states are included this additional line in the middle of 
the resonance vanishes. This effect is due to uncoupling 
of a certain phonon mode which is strongly coupled in 
the absence of the 3-phonon states. A more detailed dis- 
cussion based on a quasi-degenerate perturbation theory 
has been presented in our previous paper 

For a given distribution of QD characteristic energies 
Tiuj\^\ the absorption spectrum of an inhomogeneously 
broadened QD ensemble can be calculated. The distri- 
bution of QD characteristic energies is assumed to be 
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described by a Gaussian function 



fe,a{,huJo) 



1 



e 2<T 



where e — TiUq and a are the mean transition energy and 
its standard deviation, respectively. 

The absorption spectra of an inhomogeneously broad- 
ened QD ensemble calculated using models with different 
number of phonons are presented on Fig.[2l Again, in the 
second-order resonance area (B), only the difference be- 
tween models with up to 2 and 3 phonons is particularly 
large. After including 3-phonon states the position of the 
resonance shifts towards lower energies and the resonant 
feature (B) consists of a sharp peak only. Influence of the 
additional 4-phonon states on the second-order resonance 
area is negligible, but it is quite large if one considers the 
phonon replica of the second-order resonance, that is, res- 
onance between 1-phonon and 3-phonon states (feature 
D). The presence of 4-phonon states has similar effect 
on this area as including the 3-phonon states had on the 
area of the second-order resonance. 
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FIG. 3: Influence of QD height to width ratio on the position 
of the resonant feature (a) and on the width of the resonance 
(b). 

The influence of the QD shape (height to width ra- 
tio A) on the system spectrum in the resonance area is 
shown in Fig. [31 In a single-dot polaron spectrum it shifts 



the resonance and changes its width. For QDs with a 
larger value of A, the electron-LO-phonon coupling con- 
stants have smaller values, thus the width of the reso- 
nance decreases. As it turns out, changing this shape 
parameter does not produce strong shifts in the position 
of the second-order resonant feature. Since the position 
of the resonant feature is a nearly linear function of the 
QD height to width ratio the absorption spectrum of an 
ensemble of QDs differing not only in sizes but also in 
heights can be characterized by average of A parameter. 



IV. CONCLUSIONS 

In this paper, we have theoretically modeled the 
second-order resonance between p-shell electron states 
confined in a quantum dot and s-shell state with 2 LO 
phonons. We have studied the effect of restricting the 
computational basis on the absorption spectrum of an in- 
homogeneously broadened QD ensemble. The presented 
data confirm that including 3-phonon states is necessary 
in modeling the second-order resonant polarons. When 
considering higher order effects a theoretical model has 
to account also for 4-phonon states. 

We have also studied the dependence of the polaron 
spectrum on the QD geometry. We have shown that both 
the position of the resonance and its width depend on the 
QD height to width ratio. Due to stronger interaction be- 
tween electrons and phonons in the case of flat dots, the 
width of the resonance increases with decreasing height 
to width ratio. The character of this dependence as well 
as the dependence of the position of the resonance on the 
QD height to width ratio is nearly linear. 

The presented data confirm that the resonant feature 
is clearly visible in the absorption spectra and the system 
can be modeled with standard material parameters and 
computational techniques, which provides a way towards 
modeling of even more complicated electron-phonon sys- 
tems. 
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